Unconventional magnetization processes and thermal runaway in spin-ice Dy 2 Ti 2 7 
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We investigate the non-equilibrium behavior of the spin-ice material Dy2Ti207 by studying its 
magnetization as a function of the rate at which an external field is swept. At temperatures below 
the enigmatic "freezing" temperature T cqu ii m 600 mK, we find that even the slowest sweeps fail to 
yield the equilibrium magnetization curve and instead give a smooth, initially much flatter curve. 
For higher sweep rates, the magnetization develops sharp steps accompanied by similarly sharp 
peaks in the temperature of the sample. We ascribe the former behavior to the energy barriers 
encountered in the magnetization process, which proceeds via flipping of spins on filaments traced 
out by the field-driven motion of the gapped, long-range interacting magnetic monopole excitations. 
In contrast, the peaks in temperature result from the released Zeeman energy not being carried 
away efficiently into the bath, with the resulting heating triggering a chain reaction. 



New phases go along with new excitations, which in 
turn manifest themselves as new phenomena in experi- 
ments. In the case of the magnetic pyrochlore compounds 
RE 2 TM 2 7 , where RE = Dy or Ho and TM = Ti or 
Sn, the new phase, named spin ice, is peculiar in that 
it is highly degenerate, with a non-vanishing "zero-point 
entropy" on account of the highly frustrated magnetic 
couplings [l[ and an effective low-energy description pro- 
vided by a classical emergent gauge field |2|. 

The concomitant new excitations are - unusually for 
a three-dimensional Ising magnet - pointlike topological 
defects. These are charged not only under the emergent 
gauge field but are also sources of a conventional mag- 
netic field, whence their name magnetic monopoles [3|. 

The experiments we report and analyse here are re- 
lated to the magnetization processes under an applied 
magnetic field. This is a fundamentally non-equilibrium 
phenomenon that has a long and important history in 
condensed matter physics. It is associated with concepts 
such as the Barkhausen noise 0] . It is also related to the 
physics of "magnetic deflagration" in the case of simple 
paramagnets [5|, with links to the physics of explosives, 
for which detailed analytic theories were proposed Q ■ 

In spin ice, we find that there are three fundamentally 
different types of behavior for the magnetization process. 
Besides the equilibrium magnetization curve for infinites- 
imally (at low temperature, T, unattainably) slow sweep 
rates, we find smooth behavior for slow sweeps, which 
gives way to magnetization jumps accompanied by ther- 
mal runaway at faster sweeps. 

This richness of behavior is a consequence of the pe- 
culiar nature of the spin ice state, which requires the 
motion of thermally activated monopoles to change the 
magnetization, leading to an exponentially slow response 
rate at low T. The slow-sweep response is special in 



that it reflects the nucleation of the equilibrium (magne- 
tized) phase as in the case of a supercooled gas, with the 
crucial distinction that the nucleation happens through 
one-dimensional filaments (the paths of monopoles), so 
that the process can occur gradually even in the ab- 
sence of quenched disorder. The relevant energy bar- 
riers arise from the competition of magnetic Coulomb 
and Zeeman energies involved in creating and separating 
the monopoles. In contrast, thermal runaway involves 
physics extrinsic to the spin system, namely a "phonon 
bottleneck" , the inability of the heat bath to absorb effi- 
ciently the Zeeman energy released by the spin flips. 

In the remainder of this paper, we describe and model 
these phenomena in detail. Specifically, we study the 
low temperature magnetization process in response to an 
external [111] field which takes place between two states 
within the ice-rule manifold: the zero field spin ice state 
(SI) and kagome ice (KI) 7], defined below. 

The pyrochlore lattice consists of a cubic array of 
corner-sharing tetrahedra. In the spin ice systems, the 
spins - sitting at the corners of the tetrahedra - are con- 
strained to point into the center of one of the two tetra- 
hedra they belong to. The spin configurations in the low 
energy ice-rule manifold are those in which only two of 
the four spins of every tetrahedron point towards its cen- 
ter. Positive (negative) magnetic monopoles correspond 
to tetrahedra with three (one) spins pointing in. 

To understand the transition between the zero field SI 
and KI states, it is convenient to view the pyrochlore lat- 
tice as an alternating stack of two-dimensional kagome 
and triangular lattice planes, perpendicular to the [111] 
direction. It is possible to polarise all the spins in the 
triangular layers (more strongly coupled to the magnetic 
field on account of their local easy axes) while main- 
taining the ice rules and preserving a reduced but non- 
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FIG. 1: magnetization vs field along the [111] axis, mea- 
sured at low temperature after zero field cooling, at sweep 
rate 0.1 T/min. While above T equ ii ~ 600 mK the spin sys- 
tem reaches equilibrium within the experimental measuring 
time (as shown by a reasonable Hi/T scaling for low Hi/T, 
see inset), the lowest T curves exhibit quite unusual features. 
Notably, plateaux are followed by sharp jumps in M. Elimi- 
nating the effect of demagnetizing fields, DM, gives a crucial 
piece of information (inset): the field triggers an event that 
persists even when the actual internal field Hi decreases below 
the triggering value (the negative slope seen here is a conse- 
quence of the decrease in Hi as M increases). The details of 
the curves are remarkably reproducible. 



vanishing zero-point entropy 7]; this is the KI state. To 
further increase the magnetization along the [111] di- 
rection requires breaking the ice-rules, with an energy 
penalty of A s « 5.6 Kelvin per spin flip. As a re- 
sult, Dy2Ti2C>7 exhibits a low-temperature magnetiza- 
tion plateau at the value Mki ~ 3.3 /us/Dy for fields up 
to about 1 Tesla 0. On the other hand, at small fields 
the system's response is similar to that of an isotropic 
paramagnet and one expects a linear scaling of the mag- 
netization with H/T [§]. 

Single crystals of Dy 2 Ti 2 07 were grown by the floating 
zone method in St Andrews. The samples were cut into 
prisms of approx. 2 x 0.7 x 0.5 mm 3 with the long axis 
along [111] to minimize demagnetizing effects. The mag- 
netization was measured using a purpose-built plastic 
Faraday force magnetometer. The sample was mounted 
flat on a sapphire plate thermally anchored to the mix- 
ing chamber of a dilution refrigerator. An additional cali- 
brated thermometer was attached directly to the sample. 
The magnetometer was calibrated against a SQUID. 

Results - Fig. [1] shows the measured magnetization 
along [111] after zero field cooling, as a function of field 
(with a sweep rate v = 0.1 T/m), for different tempera- 
tures. At relatively high temperatures (above T equ ii) the 
magnetization rises linearly from zero as expected and 
reaches a plateau of approximately 3.3 /xs/Dy. As the 
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FIG. 2: M vs Hi — H — DM at 200 mK, measured for 
different field sweep rates v. Jumps only happen above 
v — 0.025 T/min. At slower sweep rates, the data seem to fall 
onto a single limiting curve. Jumps also occur in field cool- 
ing experiments; the inset shows M vs Hi at 0.1 T/min and 
200mK for different starting field cooled states (as indicated 
in each curve), together with the slow curve (at 0.01 T/min) 
and the 700 mK curve. The field-cooled jumps are contained 
within these two limiting curves, as explained in the text. 



temperature is lowered to 500 mK, contrary to the ex- 
pected H/T scaling, the magnetization rises at a slower 
rate. Below 400 mK the magnetization at first remains 
tiny until a finite field is reached, whereupon it exhibits 
a steep jump which can hardly be resolved by our instru- 
ment, recording 1 point every ~ 2 s. Further plateaux at 
intermediate values of the magnetization and subsequent 
jumps are seen before reaching Mki. With decreasing 
temperatures, both the intermediate plateaux and the 
jumps become sharper and the value of the magnetiza- 
tion at the plateaux changes (from around 2.4 /Ltg/Dy to 
2.8 /j,s/Dy). The latter is a strong indication that the 
plateaux cannot be explained by the formation of an equi- 
librium intermediate magnetization pattern, in contrast 
to the case of magnetoelastic coupling [ljj ■ 

Further information can be obtained by plotting the 
magnetization as a function of Hi/T, where Hi = H — 
DM, D the calculated demagnetizing factor based on the 
geometry of the samples (inset of Fig. [1]). As expected, 
the Hi/T scaling works only for the high temperature 
data, where it also agrees with MC simulations of the 
dipolar spin ice model (not shown). At lower tempera- 
tures the curves lie below the equilibrium curve. What is 
most striking is the fact that the jumps, plotted against 
the actual internal field, acquire a negative slope. This 
shows that they correspond to triggered events: once the 
process has started, it does not stop even though the 
internal field falls below the triggering value. This rep- 
resents a magnetic Zeeman equivalent of deflagation of a 
combustible material set off by a spark. 

Jumps in the magnetization have been observed indi- 
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rectly in neutron studies of Dy2Ti207 llj, [l2J and were 
absent from the magnetization measurements reported in 
Ref. H. Fig. [2] reconciles this apparent inconsistency as a 
consequence of the different speeds at which the magnetic 
field was swept in the two experiments. The figure shows 
the magnetization curves at 200 mK for sweep rates span- 
ning a factor of about thirty, each step progressing by ap- 
proximately a factor of two from the previous rate. They 
can clearly be separated into two groups: those showing 
one or more jumps, with sweep rates v greater or equal 
to 0.025 T/min, and those showing a continuous growth 
of the magnetization, from 0.01 T/min to 0.003 T/min. 
The changes within each group are very small, despite v 
changing by a factor four: the heights of the jumps are 
unaltered and so is the slope of the continuous curve. 

The field at which m jumps can be tuned by prepar- 
ing the sample into different initial states of magnetiza- 
tion. Shown in the inset of Fig. [5] are M vs Hi curves at 
0.1 T/min measured after cooling the sample to 200mK 
under different applied magnetic fields; a slow curve 
(O.OlT/min) at the same temperature and the 700 mK 
curve are also included. We can see that the magneti- 
zation value for each field cooled curve hardly changes 
until the curve crosses the slow rate trace. Shortly after, 
the magnetization suddenly jumps up to its equilibrium 
value at T ~ 700 mK. The slow rate curve can thus be 
considered as an approximate limit for the stability of 
the out-of-equilibrium magnetization. 
Low-temperature behavior - Focusing on field values be- 
low 0.1 Tesla, we see that the magnetization at low tem- 
perature rises slowly, reaching higher values the smaller 
the sweep rate. Qualitatively, this behavior can be un- 
derstood as a consequence of the kinetics of the mag- 
netization process in spin ice, which proceeds via field- 
driven motion of the magnetic monopoles. As their den- 
sity and speed of motion is not significantly affected 
by weak fields, the maximal rate at which the system 
can respond via this mechanism is independent of the 
value of the applied field: (dM/dt) max ss p m 10/is/ms, 
where p m denotes the density of thermally activated un- 
bound monopoles, strongly suppressed at low temper- 
ature. Here we assume a typical single spin flip time 
scale of 1 ms [13J, [l4|. Even the slowest experimental 



sweeps require processes where M changes much faster 
than (dM/dt) max to maintain equilibrium for T <C T equ ii. 
As a result, the system enters a strongly out of equilib- 
rium regime where the magnetization remains very small 
despite the presence of an applied magnetic field. 

When the sweep rate is sufficiently low, this regime 
is followed by a smooth, seemingly rate-independent in- 
crease in magnetization up to Mri (Fig. [5]). Let us 
thus consider how monopoles, the agents of magnetiza- 
tion changes, are created out of the SI ground state in 
presence of a Zeeman field. On one hand, the bare cost 
of a spin flip, which creates a neighbouring monopole- 
antimonopole pair, A s , leads to a tiny density at low T, 



say 100 mK. On the other, a single monopole can change 
the magnetization by effecting O(L) spin flips as it is 
swept to the sample surface, L being the linear size of 
the system along the field direction. 

We are thus dealing with a nucleation process where 
a spin flip out of the ground state is the part of an ac- 
tivated process akin to those giving rise to an interface 
tension between two phases in a first-order transition. 
The Coulomb attraction experienced by the monopoles 
as they move apart will then be reflected in the kinetics 
of the nucleation process. 

We emphasize that unlike a supercooled gas, where a 
nucleated droplet grows to a finite volume fraction 0(L 3 ) 
- a spin flip in SI only nucleates a filament of size O(L), as 
mentioned above. In such nucleation processes, the limits 
of long times and large system sizes do not commute but 
it is easy to see that there will be a parameter range 
where the area density of filament creation is low enough 
for their Zeeman energy release to be absorbed by the 
rest of the system - this is our slow-sweep regime [15j . As 
we will argue below, the fast-sweep regime corresponds 
to the situation where this is no longer the case, and 
thermal runaway is induced. 

The smooth magnetization curve thus provides infor- 
mation about a sequence of energy barriers, which in 
principle depend on the field direction. In particular, its 
derivative dM/dH at low temperatures is tantamount 
to a histogram of the distribution of such barriers. A 
plausible origin for these barriers can be found in the 
Coulomb energy needed to separate a thermally excited 
monopole pair. Indeed, analytical considerations based 
on the effective monopole description of spin ice (as well 
as preliminary numerical simulations) yield energy scales 
that are in broad agreement with the experimental re- 
sults for dM/dH (not shown). However, a detailed mi- 
croscopic modelling of this phenomenon is still needed 
to confirm the origin of the magnetization curve at low 
sweep rates. The [111] field direction is particularly com- 
plicated as monopoles sweeping through the system can 
become stuck when their path encounters a spin in a tri- 
angular plane already pointing along the field. Such stuck 
monopoles give rise to local stray fields that broaden the 
distribution of energy barriers around them. Moreover, 
spin flips in triangular and kagome planes have different 
Zeeman energies in a [111] field. 

Thermal runaway: the fast sweep regime - We now turn 
to the magnetization jumps, the defining distinction be- 
tween the slow- and fast-sweep regimes. We believe that 
the increased rate at which Zeeman energy is dumped 
into the system as the field sweep rate increases overtaxes 
the ability of the lattice (the phonons) to equilibrate the 
system with the bath. As a result, the sample heats up 
locally, leading to the creation of more (and more easily 
unbound) monopoles. These in turn dump more energy 
as they move in the field direction and thermal runaway 
is 'ignited' above a critical sweep rate. 
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FIG. 3: The temperature of the sample (panel a) follows 
closely the Zeeman energy release W z due to the motion of 
monopoles in a magnetic field; the latter is proportional to 
Hi(t) dM / dt (panel b). The temperature reached at each of 
the thermal runaway episodes is approximately independent 
of the sample initial temperature (panel c). 



tally through AC-susceptibility, show an approximately 
exponential slowdown below T equ ii 
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In conclusion, we have found that the magnetization 
processes in spin ice exhibit a variety of novel out-of- 
equilibrium phenomena: nucleation of the equilibrium 
phase in one-dimensional filaments; thermal runaway fol- 
lowing 'supercooling'; a distribution of energy barriers 
due to long-range interactions and a 'speed-limit' on 
macroscopic equilibrium rearrangements resulting from a 
paucity of agents to effect such changes. These are exhib- 
ited robustly (i.e. without special preparation or a careful 
choice of parameters) . Clearly, further studies of equilib- 
rium and non-equilibrium dynamics and transport - in 
particular, thermal transport - are highly desirable[l6j|. 
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To verify this, we attached an additional thermometer 
directly to the sample to measure its temperature and 
magnetization simultaneously. In Fig. EJa) and (b), we 
show this temperature for four different magnetic sweep 
rates as a function of field, along with W z cx Hi(t) dM/dt, 
the rate of Zeeman energy release. The sample is strongly 
coupled to the bath (the dilution refrigerator mixing 
chamber) and therefore the local temperature Tj should 
reflect W z . Indeed, W z and Ti trace each other closely. 
As seen in Fig. [21 for slow sweep rates the energy is dis- 
sipated gradually during the whole magnetization pro- 
cess. Once the critical sweep rate is reached (about 
0.025 T/min in this case), the magnetization jumps 
abruptly, accompanied by a strong spike in the temper- 
ature. 

This picture is further supported by the observation 
that the maximal temperature reached by thermal run- 
away is approximately independent of the initial temper- 
ature over a wide range, pointing at an intrinsic feature 
of the thermal coupling in the spin-lattice system. This is 
displayed in Fig. [31c): at T oqu ii w 600 mK, thermal equi- 
libration is efficient, resulting in a flat temperature trace. 
By contrast, for T < 500 mK temperature spikes appear 
which, crucially, never surpass T equ ii, where efficient ther- 
mal contact is clearly reestablished. This also explains 
why the magnetization always jumps up to roughly its 
equilibrium value at Tk T equ ii (inset of Fig. [2]). 

Precisely what it is that happens microscopically 
around T equ ii is an important open question, which is 
all the more intriguing because another phenomenon is 
observed there: spin autocorrelations, as probed in nu- 
merical simulations without phonons [14] and experimen- 
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